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Q, ' Abstract 

■^^ ■ We consider a parabolic version of the mass transport problem, and 

show that it converges to a solution of the original mass transport prob- 
lem under suitable conditions on the cost function, and initial and target 
^ ' domains. 



1 Introduction 



We are concerned in this paper with solutions to the optimal transport problem, 

which reads as follows. 
ly^ . Given two domains fl and 51*, and two probability measures fi and v defined 

^^ ' on them, along with a real valued cost function c defined on f2 x 57*, we wish 

r~^ . to find a measurable mapping T : f2 ^^ 57* satisfying T^^ = v (defined by 

O ■ T^KE) = ^J.{{T)-'^{E)) for all measurable E C VL*) such that 

o ' 



c{x,T{x))dfi = max / c{x,S{x))d^. 



/^ ' Under mild conditions on c and the measures /j, and u, it is known that the 



solution to this problem exists. Additionally, if 57 and 57* are subsets of i?", 
and fx and v are absolutely continuous with respect to Lebesgue measure, T can 
be determined from a scalar valued potential function satisfying the following 
equation in an appropriately weak sense: 

J det (V^ u{x) - A{x, Vu{x))) = B{x, Vu{x)), x G 57 
|t(57) = 17* 

where ^ is a matrix valued function and _B is a scalar valued function defined 
from c and the densities of the two measures. 

Under certain conditions on the domains, cost, and measures, the interior 
regularity of the potential u has been shown by Ma, Trudinger, and Wang in 



[S] , and global regularity by Trudinger and Wang in the subsequent [5] . In [S] , 
Schniirer and Smoczyk analyze a parabolic flow which is close to the correspond- 
ing optimal transport problem with the cost c{x, y) = \x — j/p. 

In this paper, we are concerned with a parabolic flow leading to the solu- 
tion of an optimal transport problem, with cost functions other than the case 
c{x,y) = \x — 2/p. More specifically, we look at solutions to the equation 

(uix,t) =logdet{V^u{x,t) ~ A{x,Vu{x,t))) 

I ~logB{x,Wu(x,t)), xen 

' G{x,\7u{x,t)) =0, x<Edn 

Mt=Q = "0 

with A and B as above, and an appropriate boundary condition G. It turns 
out that the conditions given in [9] along with a few restrictions on the initial 
condition are sufficient to guarantee long time existence to this parabolic flow, 
and convergence to the solution of the optimal transport problem as t ^ od. 
We also provide here a reference table for the notation used in this paper. 
Notation Definition Location 



u* 
h 

h* 

Y{x,p) 

X{q,y) 

A{x,p) 

T{x,t) 

Bix,p) 

B{x,p) 

Gix,p) 

Wij{x,t) 

I3ix,t) 
M 



Outer unit normal to dfl 
Outer unit normal to dfl* 
Defining function for H. 
Defining function for H.* 
Inverse of the map y i— > Vxc{x, y) 
Inverse of the map x i— > Vyc(x, y) 

Vlc{x,y)\y=Y(x,p) 
Y{x,Vu{x,t)) 

|detV^,c|- ,/,^"\, 

I x,y I g(Y(x,p)) 

log B{x,p) 

h*{x,Y{x,p)) 

Uij{x,t) - Aij{x,\7u{x,t)) 



P 

q 



h*i {yY''^ (x, y)\y=Y(xyu{x,t)) 



Gp(a;,p)|p^v«(2:,t) 

rJnBe(xo) 

sup^gaj^|V2u(a;,t)| 
^ sup^gn|V2u(x,t)| 

M^ s\xp^^QQ\wij{x,t)\ 

M^ sWx&n\'Wij{x,t)\ 

r^ {xGn\dist{x,dn) <e} 

T*, {y G n*\dist{y,dn*) <_e*} 

h Function defined from h with a c-convexity property 

h* Function defined from h* with a c-convexity property 

G{x,p) Function uniformly convex in p 

The author would like to thank Micah Warren for suggesting this 
and many fruitful discussions regarding this paper. 
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2 Preliminaries of Optimal Transport 

In this section, we recall some basic facts and definitions regarding the optimal 
transport problem, along with the key conditions from [9". 

Let H. and H,* be open, smooth, bounded domains in M", with / and g smooth 
functions on il and fl* respectively. We assume the mass balance condition. 



f= 9 (2.1) 

o Jsi* 

along with the bound 

0<A</,g<A<oo (2.2) 

for some constants A and A. We will assume the following conditions on c: 

c e C"*+"(n X TF) for some a e (0, 1]. (AO) 

We assume the mappings y i— S- ^xc{x,y) for each x G fJ and x i— > Vyc{x,y) 
for each y G H.* are injective. For any p e Vxc{x,il*) and x G fl, (resp. 
q € \'yc{n,y) and y £ Q,*) we write Y{x,p) (resp. X{q,y)) for the unique 
element of il* (resp. il) such that 



^xCix,y)\y=Y{x.p) =P 
^yC{x,y)\x=X{q,y) =9- 

We also assume a nondegeneracy condition on the cost c: 



(Al) 



det\7iyc{x,y) ^0,\fxen,y en*. (A2) 

Finally, writing A(x,p) — V'^c{x,y)\y^Y{x,p) we assume 

Dp.p^Aki{x,p)^^^.jr]kr]i > 0,Vx e r2,p e Vxc{x,Q,*),(, ± i]. (A3w) 

Additionally, we need to make the following assumptions on the domains fl 
and fl*. We write &,- hi = it—it— ■ ■ ■ tt-tt- • • • c and indicate the inverse of 
a matrix by raising its indices. 

Definition 2.1. We say that n is c-convex with respect to n* if the set Vyc(r2, y) 

is a convex set for any y G fi*. Likewise, il* is c* -convex with respect to fl if 

the set V xc{x, $7*) is a convex set for any x € fl. 

We say that n is uniformly c-convex with respect to fl* if it is c-convex and 

satisfies 

[D^i^jix) - c^'^ Cij^i{x,y)vk{x)]TiTj{x) > 5i, Vx e dfl, y £ n* (2.3) 

for some di > 0, where r is any unit tangent vector to dfl, and v is the outer 
unit normal to dVi. Likewise, J7* is uniformly c* -convex with respect to Vi if it 
is c* -convex and satisfies 

[D,y*{y) - c''''ci^,,{x,y)vl{y)]T:T*{y) > 5*, Vy £ 9^!*, xeQ 

for some (5* > 0, where r* is any unit tangent vector to dfl* , and v* is the outer 
unit normal to dfl* . 



Remark 2.2. For some fixed y G fi*, given any two points pi = Vyc{xi,y) 
and p2 — Vyc(x2,y) for xi, X2 G il, we define the c-segment with respect to 
y between xi and X2 as the inverse image of the straight line between pi and 
P2 under the map Vyc{-,y). It is clear that ft is c-convex with respect to fl* if 
and only if every c-segment with respect to any y G f2* between any two xi and 
X2 G ri remains inside fl. 

An analogous definition and remark hold for a c* -segment with respect to some 
X E fl between two points yi, J/2 G ^*- 

Definition 2.3. For a yo (z fl* and a Aq G M, we call a function of the form 
c{',yo) + ^a, 0- c-support function. 

We say that a function (j) is c-convex, if for every xi^ £ Q. there exists a 
c-support function such that 

<P{xq) = c(xo,2/o) + Ao 
(j>{x) > c(x, 2/o) + Ao, for all x ^ xq. 

We say (j) is strictly c-convex if the second inequality above is strict. 
If (j) is C^ , we say that it is locally, uniformly c-convex if 

</)„(x) - Ajj(x,y(x,V0(x))) > 

as a matrix, for every a; G 51. 

3 The Main Problem 

The symbols V, V^ , and Di will denote differentiation in the x variables, with 
specific variables indicated by subscripts when necessary. The notation u will 
indicate derivatives in the t direction. Now, by the smoothness of dfl*, we may 
extend g to all of M" so that it is C^ and satisfies the bound -I < .g < 2A. 
Writing 

T(x,t) = Y(x,Vu{x,t)) 

j(x) 
B{x,p) = \ActVl^^c{x,y)\y^Y(x..p) ■ g{Y{x,p)) 

B{x,p) ^ log B{x,p) 

G{x,p)^h*{Y{x,p)) 
Wij(x,t) —Uij(x,t) — Aij{x,Vu{x,t)). (3-1) 

where h* is a normalized defining function for 17* (meaning V/i* = v* on dfl*, 
h* = on dil*, and h* < on V,*), we consider solutions u — u{x,t) to the 
following equation: 

m(x, t) = log det (V^ u{x, t) - A{x, Vm(x, t))) 

-logB{x,Vu{x,t)), xen (3.2a) 

G(x, Vu(a;, t)) ==0, x e dQ (3.2b) 

u|t=o = -"0. (3.2c) 



We require the following conditions on uq G C^"'""(ri x {0}): 

uq is locally, uniformly c-convex (3.3a) 

< h*{Y{x, Vmo(x))) = on arj (3.3b) 

^ Tn{n)=n* (3.3c) 

where To{x) — Y{x, Vuo{x)). 

The main theorem we prove is the following. 

Theorem 3.1. Suppose that c satisfies conditions (jAOp - (jA3wP . Additionally, 
suppose that fl and fl* are uniformly c and c* -convex with respect to each other. 
If Uq G C^^"(57 X {0}) satisfies the conditions (J3.3I) . there exists a solution 
u to equation p.2p for all times t > which is C'^{Vt) in the x variables and 
C^{^+) in the t variable. 

In addition, u{-,i) converges in C'^{Vl) to a c-convex function u°°(-) as t — >■ cxd 
which satisfies the elliptic optimal transport equation: 

det{\/^u'^{x)- A{x,Vu'^{x)))=B{x,Vu°^{x)), xefl (3.4a) 
Y{-,Vu°°{-))[n]=n*. (3.4b) 

Remark 3.2. Note that if 4>{-) is differentiable at some xq G VL, and c{x, j/o) + -^o 
is a c- support function to 4> at Xi), we will have that V(j){xo) — VxcixQ^yQ). Thus, 
by the uniqueness in (jAip . we have that yo = Y{xo,'^(f>{xo)). 

4 Short-time Existence 

We will prove the existence of a solution to our equation (13.21) up to some small 
time tmax > 0. First, we will prove some auxiliary results. We follow the 
definitions for Holder spaces given in [3]. 

Definition 4.1. Write X = {x,t), Xq = (xo,to) €$7x7 for some interval I, 
and \\X\\ — max (|a;|, |t|2 ). For a function f defined on fl x I, with k a positive 
integer and < a < 1, we define the following norms and seminorms: 

ffl - V sun \D'.DlifW-fiXo))\ 
[Jjc+^cox/) - 2^ sup — —7^ , 

/r. sr . \DiDi{f{xo,t)-f{X,))\ 
(/)c"=+°(ax/) = }_^ sup — — -^ , 

Iflc+^inxl) = [f]k+a + {f)k+a 

WfWc^^in.i) = E sup|7^f^^/| + 1/1,.+,,. 

|/3|+2i<fc 

We write f G C^+^{n x /) if l|/l|c-+»(ax/) < oo- 



Remark 4.2. By looking at sequences of points in ^x I approaching d^ x I or 
ri X {0}, it is easy to see that if f G C'^"'"", the x derivatives up to order k and the 
t derivatives up to order ^ of f exist and are continuous up to {dil,x I)\j{Qx {0}) . 

Definition 4.3. 

K. = {(p £ C {n X I) I (j){-,t) is locally, uniformly c- convex for each t G /} . 
For this section, let us write 

F{x,p,r) =logdet(ry - A{x,p)) ~\ogB{x,p) 

well defined for {x,p, r) in some subset of il x M" x R"' (in particular, F is well 
defined when r = V^ <l)[x, t), p — S/(j){x, t) for e /C.) 

By Theorem 17. II below, the boundary condition p.2bl) is a nonlinear oblique 
condition. By a modification of the argument in [2l Theorem 2.5.7], we can 
show short-time existence for our equation (|3.2[) . The main difference is that 
the nonlinear boundary condition makes the set of admissible solutions not a 
Banach space, but a Banach manifold. 



Theorem 4.4. For some t„iax > 0, we can find a solution to (13.21) in the class 
C*'^(f2 X [0,tmax))- This notation means the solution u is C* in the x variables, 
and C^ in the t variable. 

Proof. Consider u, the solution to the problem 

Am — u = Auo — F[x, Vuq, V^ uq), x £ fl 
G{x,Vu{x,t)) = 0, x&dn 

^u\t=o = uo- 

By Theorem 17. 11 we have {Gp, v) > 0, thus on [0,e) for some small e, we find 

Also, by assumption G{x, \7uo{x)) = for a: G dfl, so by [H Theorems 8.8 and 
8.9] there exists a solution u G C'^'^"{il x [0,s)) (for a possibly smaller e > 0). 
In particular, ||u(-,i) — uq{-)\\(j2 < C\t\" for each t. Thus, since uq is locally, 
uniformly c-convex by assumption, (by making e smaller if necessary) we can 
ensure that m G /C. Writing w*-' for the matrix inverse of Uij — Aij{x,'S/u) we 
see that w*-' is positive definite for i < e and w^-'^i^j > C^\£,\'^ for some C^ > 0. 
Define 

/(x, t) := F(x, Vu, V^ u) - u{x, t). 

Note that / G C"(r2 x [0,e)) and f\t=o = 0. 

We now define the following sets for < (5 < |: 

i3,-:={<^GC2+f(f7x[0,£)) I U~u\\c^^f(n.^^^,))<^} 



and 

Bl :=C7(0x [0,e)) xS 

where 

Clearly, Bg is a C^ Banach manifold with charts diffeoniorphic to subsets of 
C72+f(fix [0,e)). 

We also claim that S is a Banach manifold. Consider the map H : C^+? (il x 
{0}) X Ci+7(af7 X [0,e)) ^ C^+T{dn X {0}) given by 

H{v, w) := G{x, \7v{x)) ~ wix, 0). 

We see that the differential 

D{vo,wo)H{v,w) = {Gp{x,Vvo{x)),\7v) - ■w{x,0) 

is clearly onto (7^+^(90 x {0}) for each (uo,wo) G B. Also, we consider the 

map P : C^+f (O x {0}) x C^+f {dQ x [0, e)) -^ J\f{D^^^^.^„^H{v, w)), the kernel 
of -D(^o,^o)iJ(w,w), defined by 

P{v,w) := (w,w + (/)„_„,) 

where 

</>«,«) (a;, t) = {Gp{x,Vvq{x)),'S/v{x,0)) —w{x,0), for ah {x,t) e dVl x (0,£). 

Since P is a continuous linear projection, M{Diy^y^^\H{v,w)) splits the full 
tangent space of C^+7(f7 x {0}) x C^+^{dQ. x [0,e)) at (uo,wo) and we have 
that {vq^wq) is a regular point of H. Thus, B is the inverse image of a regular 
value of H , so by [lOl Theorem 73C] is also a Banach manifold with tangent 
space at {vo,wq) equal to A/'(I?(.u|-,,^g)i/(f, w)). Hence Bg is also a C^ Banach 
manifold. 

We now define the map 

^■.Bs^ Bl 

by 

$(u) = (F(a;, Vu, V^u) -u, u\t=o, G[x,Wu)). 

By the uniform continuity of Qj, we can see that for a function v G G'^{Vl), if 
V'^ V — A{x, Y{x, Vu)) is positive definite, it remains so on a small neighborhood 
of V in the G'^{fl) norm. Hence, if S is sufficiently small, any convex combination 
of u and (j) ^ Bg will remain in /C. Thus we find that the differential of $ at m 
is given by 

£ii$(^) = (w'^'((/)y - Dp^Aijix, Vu)(i>k) - Dp^B{x, Vu)4>k - 4>, 
</'|t=o, {Gp{x,Vu),V4>)). 



Since by above, the tangent space to i3| at <J>(u) is C^(r2x[0,e))xA/'(D($(ii))i?), 
we can use O Theorem 5.18] to find that D^^ is bijective on its tangent space 
(the condition that a pair be in the kernel of the differential of H is exactly the 
required compatibility condition between the initial and boundary conditions 
to obtain short time existence for a linear parabolic equation). Hence, by the 
inverse function theorem for Banach manifolds, $ is invertible on some small 
neighborhood of u G Bs- 

Now, taking a smooth cutoff function rj in the variable t such that 

Jo, 0<t<5 

with \?)s\ < y, consider the function fs{x,t) := ris{t)f{x,t). Since f\t=o = 0, 
we can make sup \fs — f\ small by taking 6 small. On the other hand, we 

XeOx[0,e) 

calculate that 

\js - /lcT(ox[o,£)) - I'' let (ox[o,<5)) + ^\vs - ^\cf {nx[s,2S)yj\c^ {nx[s,2S)) 
where C is a constant depending only on fl and ||/||c°(r2x[o,e))- Since 

|/>,i)| = |/(x,t)-/(x,0)|<CP, 
we have 

|?| ., |/(x,^)-/(.To,to))| , \f{x,t) - f{^oM))\ 

^c f \f(^'^) - fixo^to))\ tu .-, fi . ^^|3V 

< sup — |/(a;,i) -/(a;o,io))| 

\ |X Xq I / 

]/(M)_i/(^o^io))|| 1^ 

+ &UP — ^ \t-tQ\ 8 

while by a similar calculation, 

1^^ - llcf (nx[,..))l/lc*(ax[,..)) < C (-P^^^ff^^) (-5^ +^-) 

<cfsupM-:^)(,f+,f) 

Thus for a small enough 5 > Q, rjsf \s arbitrarily close to / in C^{Q. x [0,e)). 
This gives a u G C^+^(il x [0,e)) such that 

F{x,Vu,V^u) ^rjsf, xeQ 
Gix,Vu{x,t)) = 0, xedn 

u\t^o = Uq. 



Hence, by taking tmax < ^, we see that u is the solution to our original prob- 
lem pT2)) on [0,tmax) in the space C^+7(ri x [0,t„iax))- Thus by Remark (|4?2|) . 
we have the desired regularity. D 

From here on, we assume that u exists with this regularity for t G [0,tmax)- 

5 Preliminary results 

We will show a few preliminary results before proceeding to show the estimates 
necessary for long term existence. 

Remark 5.1. By implicit differentiation, we see that DiT^ {x, t) = c-'' Wki{x, T{x, t)) 
and thus det DT = -f^^e" from (I3.2ap . In particular, det DT ^ for all x G Q 
and t £ [0,tmax)- We will use this fact frequently. 

Remark 5.2. We will also make use of the following formulas obtained by a 
simple differentiation: 

C^ — Cip.qC C 

J^'^ - ^c . J,Pf,q,k 

Lemma 5.3. If uq satisfies the conditions (J3.3I) . the solution u{-,t) to (13.21) 
is both locally, uniformly c-convex and strictly c- convex for < t < tmax- In 
particular, Wij remains positive definite for < t < tmax ■ 

Proof. By the assumption on the initial condition uo, we see that Wij remains 
positive definite for at least some small time. Then, since det Wij — Be^, by 
conditions (|A2[) and (|2.2p wc see that Wij cannot have for an eigenvalue, i.e. 
it will remain positive definite as long as the solution exists, proving the claim 
of local, uniform c-convcxity. 

Now, suppose that for some t, u{-,t) is c-convex but not strictly c-convex, 
i.e. there exists xi ^ X2 such that c(-, y) + A, for some y and some A G R is a c- 
support function to u(-,i) at both xi anda;2. Define wi(a;,t) = u{x,t) — {c{x,y) + 
A). By local, uniform c-convexity we see that 'V'^ vi(xi) = {wij{xi)) > 0, while 
by the definition of T and Remark 13.21 we have that Vwi = at x = xi. This 
implies that vi{xi,t) = is a strict local minimum on some neighborhood of 
xi. By the continuity of the map X{q, y) in ()A1|) . we can pick a point xq on the 
c-segment with respect to y between xi and X2 that lies in this neighborhood of 
xi. The c-convexity of fi with respect to 17* ensures that xq G fl. Since u{-,t) 
is c-convex, there is a yo G ^* and a Aq G R so that c{x, yo) + Aq is a c-support 
function at xq. In particular, 

c{x, y) + \ = u{x) > c{x, yo) -f- Aq at a; = xi, X2. 

Since we assume condition (IA3w[) . we can use [H Theorem 3.2], to show that 
Ao - A < 4>{y) < c{xo, y) - c{xo, yo), where 

(j){y) :=min{c(a;i,y) ~ c{xi,yo),c{x2,y) - c{x2,yo)} 



(recalling that xq lies on the c-segment with respect to y between xi and 12)- 
We make note that, due to the fact that the sign of the potential function in [4] 
is the negative of ours, and the x and y variables are reversed in the theorem, 
the statement becomes: 

(j){y) < c{x0,y) -c{x0,yo) 

for any xe (z ^ lying on the c-segment with respect to j/q connecting xi and X2 
and any y e fl* . But then, u{xQ,t) = c(a;o, yo) + ''^o < c{xo,y) + X, contradicting 
that xi is a strict local minimum for vi. Thus m(-, t) is actually strictly c-convex. 
Now consider the function 

v{x, xo,t) ^ u{x, t) - {c{x, T{xo,t)) + u{xo) - c{xo,T{xii, t))) 

defined on fi x fJ x [Q,tmax)- Take t < t2 < t,nax, then {wij{x,t)) > CSij as a 
matrix, for some C depending on ^2 but independent of x and t. Then, since 
Wij {x, t) = Vx v(x, XQ,t) dX X — xq, we have that 

v{x, xo, t) > v{xq, Xq, t) + {Vxv{xo, Xo, t),x ~ Xq) 

+ {Vlv{xo,xo,t){x - xo),x- xo) - \\v\\cA^ - ^o? 
> + + C|a: - xoP - hWcAx - xo\^ 
>0 

as long as < |a; — xq] < 5, for 5 depending only on C and 

||w||c3=sup sup \dl^^^^^v{xQ,x,t)\. 
^■j,k (x,t)enx[o,t2] 

Now, if |x — xol > 6, by the strict c-convexity of uq we have 

v{x,xo,0) = uo{x) - c{x,To{xo)) ~ uo{xo) + c{xo,To{xq)) 
> K{xo) > 



where K{xa) '■= inf w(a;, a;o, 0). By the continuity of uq, To, and c we have that 

X 

K(xq) is actually a continuous function of xq. Thus by boundedness of fi, we 
have that mi_K{xo) > Kq for some constant Kq > 0. Now by the continuity of 

xo£fl 

V, for some sufficiently small £ > 0, we have that \v{x, xo,t) — v{x, xq, 0)| < ^. 
Thus, for < i < e, we have 

u(x, xo,t) > 0, for < |a; — xqI < S 

v{x, Xo,t) > v{x, Xq, 0) — \v{x, XQ,t) — v{x, Xq, 0)| 

> Kn — — > 0, for la; — xnl > 5. 

- 2 2 

We can see this implies that u{-,t) remains strictly c-convex for < i < e. 
Taking ti — sup{i < imax|'u(',^2) is strictly c-convex for ^2 < t}, by the argu- 
ment above, ii > 0. If ii < tmax we have that u(-,ii) is c-convex but not 



10 



strictly c-convex, which is a contradiction. Thus we have that ii — tmax and 
u{-,t) remains strictly c-convex on < t < t„iax- 

□ 

Corollary 5.4. Under the assumptions of Lemma \ 5. S\ above, the map T{x,t) = 

Y{x,\Iu{x,t)) is one-to-one on Q for each <t < tmax- 



Proof. Suppose that y — T{xi,t) — T{x2,t) for some xi ^ x^. By Remark [ 
we have that c(x^y) + Ai and c{x,y) + A2 are c-support functions to u{-,t) at 
xi and X2 respectively, for some Ai, A2 G M. Thus 

c{x2, y) + Ai < u{x2,t) = c{x2,y) + A2 
c{xi,y) + A2 < u{xi,t) = c{xi,y) + Ai 

and hence Ai — A2. This contradicts the strict c-convexity of u{-,t), and we 
have that T{-,t) is a one-to-one function for each t < tmax- CH 

Corollary 5.5. Under the assumptions of Lemma \5.3[ T{fl,t) = fl* for < 

Proof. By c-convexity and the conditions ()A1|) and (|A2|) . Vl and il* are homeo- 
morphic to convex sets, and thus the unit ball in K", and also their boundaries 
are connected. Since the argument here will be entirely topological, we may com- 
pose T with the appropriate homeomorphisms, and assume Vl = Vt = _Bi(0). 
The boundary condition (j3.2bp implies that T{dVL,t) C dVl* for each t. T is 
continuous on dVl, and by Corollarv l5.4i is one-to-one there. By Remark 15. 1[ 
{T)~^ is also continuous on T{dQ.,t) and thus T{-,t) is a homeomorphism on 
9ri. Hence, 

T{dn, t) = dn* for each t. (5.1) 

Suppose there is some io < tmax such that T(r2, to) ?^ ^*- That means, for some 
Xq G Q., T{xo,to) ^ fl* . However, since T(a;o,0) G fl* by assumption, by using 
the intermediate value property on the continuous, real valued function (j){t) = 
|T(a;o,t)|, there must be a i G (0,to) such that T{xQ,t) G dil.* , contradicting 
that T{-,t) maps d^ one-to-one and onto dfl*. Thus, T{U,t) C fi* and hence 
T(n, t) cTF for any t < tmax- _ _ 

As above, T(-, t) is a homeomorphism on Q, as well, so we actually have T(yi, t) = 
H*, t < tmax- But since T{dn,t) — dfl* , and since T{-,t) is one-to-one on fi, 
we see that 

T{n,t)^n*,t<tmax. 

n 

Remark 5.6. By the above Corollary, we can now see that the solution u will 
be independent of the extension that we chose for g outside of 57* . 



Lemma 5.7. Under the conditions of Corollary \5.4\ define the c-transform of 
u{x,t) as u*{y,t) ~ c{x,y) — u{x,t), where y — T{x,t). Also, define T* {y , t) = 
X{Vu*{y,t),y). Then, 



11 



1. T*{;t)^T-\;t). 

2. u* satisfies the following flow equation 

' u*{y,t) =\ogdei{V^u*{y,t) ~ A*{Vu*{y,t),y)) 

-\ogB*{Vu*{y,t),y), yen* 

hiX{Vu*{y,t),y))^0, y & dn* 

^U*\t=Q = (uq)* 

for < t < tjnax where 

(uo)*(To(x)) = c{x,Tq{x)) - u{x) 
A*{y,q)^^lc{X{q,y),y) 

9{y) 



B*{y,q) = |detV^yC(a;,y)U=x(g,y) • 



f{X{q,y)) 



for f extended to he in C^(R") with hounds j < f < 2 A, and h a defining 
function for fl. 

Proof. Since we know that T(-, t) is invertible for each t, we may fix a i and write 
u*{y,t)^c{{T)-^{y,t),y)-u{{T)-\y,t),t)andx = iT)-^{y,t). We will drop 
all t for ease of notation here. Differentiating both sides in y and remembering 
the definition of T gives us 

Vu*iy) ^ DiiT)-') ■ [V.c((r)-i(y),y) - VuiiT)-\y))] + V yc{{T)-\y),y) 
- D{{T)-') ■ [\/.,cix,T{x)) - \/u{x)] + \7yc{{T)-\y), y) 
= \7yciiT)-\y),y). 

Thus by the uniqueness in assumption (jAfp . we have that T*{y, t) = {T)^^{y, t). 
Next, we differentiate both sides of the relationVu*(y) = \7 yc{{T)^^ {y) , y) to 
obtain 

V^ u* = V^ ciiT)-\y),y) + DiiT)-'){y) ■ Vl^ ciiT)-\y), y) 
^Vlc{T*{y),y) + [D{T)]-\^T'iy) ■^lyc{T*{y),y). 



Rearranging and taking determinants of both sides, and using Remark 15. li we 
obtain 



7lu* - V^c(T*(,),,)) = e-^"*^^»7(f§;^ ■ IVI,. 



Now, differentiating u*{y,t) — c{{T) ^{y,t),y) — u{{T) ^{y,t),t) in t, we have 
u*{y,t) = D{{T)-Hy,t)) ■ [V,c{{T)-\y,t),y) - \7u{{Ty\y,t),t)] 
^u{{T)-\y,t),t) 
= -u{T*{y,t),t). 

Combining the above, we obtain the desired equation for u*. D 
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6 Estimate of Vn 

Theorem 6.1. As long as a solution to the equation (j3.2p exists on a time inter- 
val [0,tmax), |Vu| < C for some C > depending on Q* and c, but independent 

Proof. By Corollarv l5.5l and the boundedness of 51*, \T\ is bounded independent 
of i. From the relationship ^xc{x,y)\y^T{x) = ^u{x,t) and the boundedness of 
IVj-cl the estimate is immediate. D 



7 Obliqueness of the Boundary Condition 

Let p and i/* be the outward unit normal vectors to ft and ft* respectively. 
Writing G{x,p) = h*{Y{x,p)), the boundary condition (I3.2b[) can be written 

G{x, Vm) = for a; e dQ. 



Theorem 7.1. As long as the solution to (j3.2p exists, we have obliqueness of 
the boundary condition, i.e. 

{Gp{x, Vu),u{x)) >0,xedn,t>0. (7.1) 



Proof. First we note that h*(T{x,t)) = for x G dVt, so for any r which is 



tangential to 957, we find that hIDiTWi = 0. Also, since ft,* < in 57, we find 



from this 

h*iD,T' = h*lc'^''wk^ = xy^ (7.2) 

or equivalently, hld^'^ = xu;*'^j/fe for some x > 0. Writing /3(a;,t) = VpG(a;, p)|p=VM 
Gp{x,Vu), we see that 13k — h'^^Y^Jp^yu — /i;*c'''^|p=vn, hence 

(/3,i^)=XW^''^fe^/. (7.3) 

By Remark 15.11 det DT ^ 0, and hence x > 0. By Lemma [5.31 vf^ will 
remain positive definite as long as the solution exists, hence the desired oblique- 
ness. D 

8 Estimate of u 

Theorem 8.1. We have the estimate 

min-u < u{x,t) < max-u, \/x E fl, <t < tmax- 

Proof. Fix some ti < t^ax- By differentiating p.2ap in t and writing v(x,t) = 
u{x,t), we find that 

Lv = w^^Vij - Dp^AijVk) - Dp^Bvk - u = 
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while differentiating the boundary condition (j3.2bp in t, we see that 

t;^ = on dil. 

Let ft. be a normahzed defining function for il, ie. Vh = i/ and ft, = on 90, 
and ft < on f7. Define vi{x,t) = v{x,t) — eh{x) — Cit for some fixed e > 
and a constant Ci > to be determined. By Lemma |5.3[ w^^ remains positive 
definite as long as the solution exists, so for any fixed e > 0, on fl x [0,ti] we 
have 

Lvi = Lv- e[w'^{hij - Dp^Aijhk) - Dp^Bhk] + Ci > 

for the choice of Ci = Ce{l + M) with a constant C depending only on bounds 
on |Vft|, iV^ftl, iVpAyl, and \VpB\, and with 

M= sup iv[w'\x,t)). 

(a:,t)eOx[0,ti] 

Thus, the maximum principle implies that a maximum for vi must occur on the 
parabolic boundary of J7 x [0, ti] 

Now if x e dn, we have {v - eh - Cit)p = -e(Vft, (3) = -e{l3, v) < 0, by dZH). 
However, if vi has a maximum at some x G dQ, we have that V(f — eh — C'lt) = 
XV for some x ^ 0- This implies that {v — eh — Cit)^ = xW, i^) > 0, a 
contradiction. Thus, the maximum for vi must occur when t = 0. Recalling 
that ft < on fi we have 

v{x, t) < {v{x, t) - eh{x) - Cit) + Cit = vi [x, t) + Cit 

< sup Vi {x, 0) + Cit < sup v{x, 0) + £ sup |ft| + Ce{l + M)tmax- 

By letting e — ?► we obtain the bound 

v{x, t) < supti(a;, 0) 
x£n 

for t < ti, and since ti < t^ax was arbitrary, we obtain 

sup v{x,t) < supv{x,0). 

(a;,t)Gnx[0,tmax) xeQ 

The lower bound is similar. D 



9 Uniform Obliqueness of the Boundary Condi- 
tion 

By Theorem lA.ll we can choose ft* in a way such that ft* < on F*. = 
{y e f]*| dist(j/,arj*) < £*}, ft* = and Vft* = v* on dn* , while we have 

[A,ft*(2/) - c'^^'cuj{x,y)Dkh*{y)M, > 6*\^\' (9.1) 

for any a; S fi, y £ F*. , and £, G M". Additionally, we take G as defined in 
Corollary ED 
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Remark 9.1. Note that the boundary condition (j3.2b[) 

G(x, Vii) ^Q ondVL 

is equivalent to 

G(x,Vu) =0 ondn. 

Also, since \7h* — V/i* — v* on dfl* , we have that Gp{x, Vm) = Gp{x, Vm) = 
fi(x,t) ondn*. 

Theorem 9.2. We have uniform strict obliqueness of the boundary condi- 
tion (|3.2bp for t e [0,tmax), i-e. 

{Gp{x,\7u{x,t)),i^{x)) > C > 0, X edil, < i < t^ax 

where C depends on Q, f2*, B, c, and uq but is independent oft and t„iax- 

Proof. We take ti < tmax, and then find a C as above that is independent of 
ii. Continuing the calculations from Theorem 1 7. 1[ we find 

WMc''''h:c'''h* = X^i{(^'''h*) = X^iPi = x(/3, ^)- (9-2) 

Since % ^^ 0, we can now combine (|7.3p and (19.21) . and use Remark (|9.ip to 
write 

= (w'''vuVi){Wrane^^'Kc^'^h*) (9.3) 

for a; G 9r2. We will proceed by bounding the two terms on the right from 
below. 

Define the linearized operator by, 

L0(a;, t) = -0(x, t) + w'^ [x, t){<j),j {x, t) - Dp,A,j {x, T{x, t))<j)k {x, t)). (9.4) 

Let (cco, io) be a point where (/3(a;, t), i'{x)) achieves its minimum on dQ x [0, ii). 
We also define F(a:,p) = {Gp{x,p),iy{x))~KG{x,p) a.ndv{x,t) = F{x,Vu{x,t)), 
where v{x) is an extension of the outward normal to a neighborhood of dD. 
using the function h. Since G{x,Vu{x,t)) — h*(T{x,t)) = for a: G dVl, v{x,t) 
restricted to dfl x [0,ii] achieves its minimum at the same point {xQ,to). Let 
^e{xo) = ^ni?e(xo) with £ chosen small enough so that h < oni}^{xo). Then 
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we calculate 

Lv = -Fp^iik + w*-' {Fxix, + 2Fx.p^Ukj + Fp^Ukij + Fp^p^UkiUij 

= Fp^Luk + w*^(i^xix, + 2Fx.p^Ukj + Fp^p^UkiUij - Dp^AijF^^) 
= Fp,{w'^Dx,A, +D.,,B + Dp.Buik) 

= Fp,{w'^ D^^Aij + Dx,i3 + Dp,i3{wik + Aik)) 

+ w^^ {F^iXj + 2Fx^p^ {wkj + Aki) + Fp^p^ (wkiWij + 2wkiAij + AkiAij) 

Pm ij Xjj-i_} 

< C(l + tr (w*^)) + 2Fx,p^w''wkj + w^^ Fp^p^WkiWij + 2Fp^p^w'^ wmAij 

+ Fp,{Dp,B)wik 

< C{1 + tr {w'')) + Fp,p,Wki + Fp,[Dp,B)wik (9.5) 

where C is a constant depending on il, i7*, k, uq, the cost c, and B. The 
expression for Luk comes from differentiation of the equation p.2ap . 
By Corollary I A. 21 we have that Gp^p^^kS,! > (^2 1^1^; so we have 

Fp^p,Wki < Dp^p^{l3,iy)wki - KGp,p,Wki 

^iCo-S;K)tT{w,j). (9.6) 



Thus by choosing k large enough, we will have by 

Lv < C(l + tr {w'^)) < Ctr {w'^). (9.7) 

The second inequality comes from the arithmetic geometric mean inequality, 
combined with the fact that 

detw'^' = (detw^j)-^ = (B + e")-i > ^ 



from bounds on /, g, V^ y c and Theorem 18. II 
By Theorem I A. 1 1 we find that 

Lh = w'^iD.jh - c''''c,j,i{x,T)Dkh) > Soti (w*^) (9.8) 

We fix a point j/q G ^* ^-nd consider the function 

e{x, t) ^ {v{x, t) - w(xo, to)) + (a, Vyc(a;, yo) - Pa) 

+ Ci\x - xo\^ -^ C2h (9.9) 

where po = Vyc(a;o, 2/o)i Ci and C2 are positive constants, and a G R" are to be 
determined. The idea is that we will adjust a to obtain 6 > on ^^(cco) x {0}, 
Ci to obtain 9 > on d{Vt^{xi))) x (0, tmax], and C2 to bound the term involving 
tr (w*-') in Lv so that LQ < everywhere. Calling the map ^(x) = \/yc{x, yo), 
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which is a difFeomorphism by conditions (|A1[) and (jA2p , we have by the uniform 
c-convexity of il that J7i = ^(O) is a strictly convex set in R". Hence, there 
is a supporting hyperplane to fli at po = ^(2^0)1 smd we may assume that 
the normal direction to the supporting hyperplane away from fJi is given by 
e„. Write vi{p,t) = w(^~^(p),t), P{p) for the projection of p onto di^i in the 
e„ direction, and Po{p) for the orthogonal projection of p onto the supporting 
hyperplane to ili. Then, since vi{p, 0) — wi(po, to) > for p G di^i we have that 

vi{p,0)-vi{po,to) > vi{P{p),0) ^ vi{po,to) - sup\\7vi{-,0)\\p ^ P{p)\ 
>-snp\Vvi{;0)\\p-Po{p)\ 
> -{a,p-po) 

if a is a sufficiently large multiple of — e„. Additionally, with this choice, {a,p — 
Po) ^ on the half space defined by the supporting hyperplane to ili at poi 
hence also on fti, and thus {a,'if{x) — po) > on n^{xo). By positivity of the 
terms —h and |a; — xo\'^, we find that > on il.g{xo) x {0} as desired, and a 
depends only on ui at t = 0. Next, on d{ile{xo)) x {0,tmax] we have that 

e > Ci\x-Xo\'^ + {v{x,t) -v{xo,to)) 
>0 

for Ci large depending on e and upper bounds on |w|, which in turn depend on 
h* and c. 

Finally, since 

L{{a,Vyc{x,yo) - pa)) = w'^'ci^^kcik - w'^^ (Dp^Aij)ci^kak < Ctr(w*-') 

and 

Li\x - xoH = tr(u;*^') - w'\Dp,A,j){x - xo)k < Ctr^^') 

we calculate from (19.71) and (19.81) that 



Le<CtT (w'^) - S0C2 tr {w'^) < 

for C2 sufficiently large. Thus noting that 9(xo,io) — 0, by the comparison 
principle we find that (VG, v) < ot 

{Vv,u) <C\{a,v)\ +C2\{yh,v)\ <C. (9.10) 

Now we calculate at {xQ,to), using the formulae from Remark (j5.2p . 

= /iL(Ar™)c'''=z/fc + /ir(cf + c';tAr'")^^fc + /^rc'-'Ai^fc 
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Multiply by ft.*c'''' and sum to obtain 

+ c"'''x^pc'''=z.fc(/^r™ - hX''cs,mi), by dlSD 
> SoJ^K^'-'l' + xSoJ2\^,c'^-'\^, by Theorem ED 

i i 

>5o^|/^;c*1'>Co, bydUl). 

i 

Since w(-, to) takes its minimum on dfl at x — .tq, Vw(a;o, to) = ao'^l^^o) for some 
ao < C from (igjU)) . Thus at (a;o,to), 

ao{P,v) + Khlc''^'h*c''^Wkr = /i*c'?''(aoi^» + k/i;c''''=u;h) 

= /^;c'?''(aoJ^, + KA(/i*(r))) 



Now if ao < 0, we can throw that term away and obtain the desired bound. 
If ao > we have two cases. If (/3, i') > |^ > |^, again we have the desired 
bound already. Otherwise, we see that 

Kh*c''''h*c'''''wM > -aa{l3,y)+Co 

>^. 
- 2 

In this last case, we continue by estimating w''^ViUj from below. To do this 
we use the transportation problem in the opposite direction. 

Taking u* as in Lemma 15.71 and again taking h and h* constructed in The- 
orem [XTl we find that for x E dfl, 

{P,y){x,t) = c''-\x,T{x,t))hi{x)hl[T{x,t)) = {fi\u*){T{x,t)), 

if we define 

where G*{q,y) = h{X{q,y)). This implies that if we take (xo,io) again as the 
point where the minimum to (/3, J^)(x,t) occurs, {13*, v) has its minimum over 
dn* X [0,ti] at iyo,to) = {T*{xo,to),to). If we write wli{y,t) = uli{y,t) - 
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c ,ki(T*{y,t),t), we find that for {x,y,t) = {x,T{x,t),t), we have c^'''w*^. = 
DjT" = (DjT')-^ = {c'''wij)-^ = w*'q,j. Hence 

wlic^'^c^^^ = w'J' (9.11) 

u-Hc'^-'^c'J = w"^' (9.12) 

and 

wlic''^'^c'^^{T*{y,t),y)vmVn{T*{y,t)) = w'^ {T*{y,t),y)v,v,{T*{y,t)). 

Thus with a similar proof, we estimate u'^,;C*'''™c''"^'m^'„ from below, and hence 
w'^^ViVj and by (|9.3p the lower bound of {P^v) is established. D 

10 Interior C^ estimates 

Theorem 10.1. For each t G [Q,tmax) we have the bound 

sup\\7^u{x,t)\<C{l + M) (10.1) 

where M — sup |V u{x^t)\, for a constant C depending on Vt, Vt* , B, c, and 

xedn 
uq, but is independent oft. 

Proof. Assume again, that ti < tmax- Redefine the linearized operator as 
L(j){x,t) = -(j){x,t) + w''\x,t){4)ij{x,t) - Dp^Aij{x,T{x,t))(j)kix,t)) 

and let 

v{x,t,() =\og{w^^{x,t)) +a|Vu(a;,i)p 

for some fixed constant a, where ^ is a unit vector. We differentiate p.2ap twice 
in X, in the ^ direction to obtain: 

w'^[D^jU^ - D^A^j - {Dp,A^j)Dku^] - u^ = D^B + {Dp^B)DkU^ 

and since D^w^^ ~ ~w'^''{D^wik)w^^ , 

w^^[DijU^^ - D^^Aij - 2{{D^Dp^Aij)DkU^) - {Dp^p^Aij)DkU^Diu^ 

- (Dp^Aij)DkU(^(} - u^^ - w^''w^''D^WijD^wik 

= D^^B + 2{D^Dp^B)DkU^ + {Dp^B)DkU^^ + {Dp^p^B)DkU^Diu^. 

(10.2) 

Now fix {x, t) and let {w*^} be a set of orthonormal eigenvectors for Wij{x, t) with 
eigenvalues A^. We write {v^)i to denote the rtli component of v^. Using (jA3wP 
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we obtain at {x,t), 
p 

q^p r^p 

p q or r—p 

p,q 

= -Ctr(wij). 

Using this and (|10.2l) we calculate 

Lu^^ = w*^(AjWcC ^ Dp^Aij)DkU^^) - ii^^ - {Dp^B)DkU^^ 

= w'^w^'^D^WijD^wik + w'^{D^^A,j + 2{{D^Dp^A,j)DkU^) 

+ {Dp,piA,j)DkU^DiU!=) 

+ D^^B + 2{D^Dp,B)DkU^ + {Dp,p,B)DkU^Diu^ 
= w'^w^'^D^WijD^wik + w'^[D^^A,j + 2{{D^Dp,Aij){wki + Ak^)) 

+ {Dp,p,A^j){wk^wi^ + 2wk^Ai^ + Ak^Ai^)] + D^^B 

+ 2{D^DpJ)iwki+Ak^) 

+ {Dp^p,B){wk^wi^ + 2wk^Ai^ + Ak^Ai^)) 

> w'-^w'^D^WijD^wik +w''^{Dp^p,A,j)wk^wi^ 

- C(l + tr {w'^) + tr {w'^)tT {w,j) + tr (wy) + tr (wijf) 

> w'-'^w^'^D^WijO^wik 

- C(l + tr {w'^) + tr K^') tr (w.^) + tr (wy) + tr (wy)^)- 

Also, 

LA^^ = w'^[D^^^^A^^ + 2{D^^p^A^^)ukj + {Dp^p,A^i^)ukiUi, - Dp^A,j{D^^A^^)\ 
+ Dp^^A^^[w''^ {ukij - {Dp^Aij)uki) - Uk] - {Dp^B)DkA^^ 
= w'^[D.^^.^^A^^ + 2{D^^p^A^(){wkj + Aki) 
+ {Dp^p^A^(){wkiWij + 2wkiAij + AkiAij) 
- Dp,A,j{D.,,A^^)] + Dp,A^^[w'W.,,A,j + D,,B 
+ {Dp,B){wik + Aik)] ~ {Dp^B)DkA^^ 
> -C(l + tr {w'^) + tr (w*^) tr (w,j) + tr (w^)) 
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where the expression for w'^ (u^ij — (Dp^ Aij)uki )—Uk is from differentiating p.2ap . 
Thus we have 

Lw^(^ > w'^^w^^ D(^WijD^wik 

- C(l + tr (w*^) + tr (w'J) tr (w^) + tr (w^) + tr {wijf). (10.3) 

Differentiating v, we have 

DiWff 
DiV = h 2auiuii 

DijWff DiWffDjWff 
DrjV = — ^-^^ , -„ + 2auijUii + 2auiuuj 



and 



Wee 
V = — ^-^ + zaukUk- 



Suppose that v takes its maximum at some point (a;o,io) in (^ x (Oj^i)) U 
{ft X {ti}) and some £,. There we have V^w < 0, Vw = 0, and v > 0, hence by 
differentiation of the equation p.2al) , writing Uij — Wij + Aij again, assuming 
w^^ > 1 and using (|10.3p we have 

0>Lv^ ^ ; ^\/ ^^ + w'^2auijuu 

+ 2aui[w''^ {uuj - {Dp^Aij)uki) ~ ui ~ {Dp^B)uik] 

> ^ — —^ , \„ + w''2awijWii + 2auiD^^B 

- C{w^()-^{1 + tr {w'^) + tr (u;*^) tr {w^j) + tr {w,j) + tr {w,jf) 

-Ca-Ctriw'i). 
For the first two terms above, we change coordinates so Wij is diagonal at (xq, to) 
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> 



with C = ei and calculate 

^■^^ w"[(Awii)^ + 2(uii, - DiAi, - M^ii + D^All){D,wn + um - DiAu ~ A^n + A^n)] 
E^>l t«"[(Awii)^ + 2(AAii - DiAu){2D,wii + AAii - DiAu)] 
J2,>i w"[(Awii)2 + 4(AAii ~ AAh)(Awii) + 2(A^ii ^ DiAuf] 
E,>iW'"[-2(AAii- AAi,)"] 



> 



{wiiy 



2 



-tr(«;'J)C(l+tr(w,,)) 
(z«ii)2 
> -CtT{w'^). 

Hence we have 

Ca + Ctr {w'^)) > (2a - C) tr (wy) 

so choosing a large enough, we obtain at (xo,to), for any e > 

Ce + e sup tr {w''^ ) > tr (w^ ) . 
nx(o,ti] 

Now again, using the solution u* defined in Lemma 15.71 using (j9.1ip , and ap- 
plying the same calculations, wc find that 

trK^') = trKic'='V''^') 
<C7tr(w*,) 

< Ce + Ce sup tr {wij). 
nx(o,ti] 

Thus combining the above, for s small enough we obtain 

sup tr (wij) < C. 
nx(o,ti] 
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If the max for v occurs at xq G dO, or to = 0, we can simply add the terms 
sup|V^ m(x, 0)1 + M to the right hand side of the estimate, and we absorb the 

first term into C by allowing its dependence on uq- CH 

11 Boundary C^ estimates 

By tangentially differentiating the boundary condition h*{T(x,t)) == on dil, 
we see that U/St = for any r tangential to dfl. 

Theorem 11.1. For each t e [0,tmax), 

upp{x,t) <C(1 + M)^ 

for some C that depends on fi, f7*, B, c, and Uq, but is independent oft, where 
M = aup\V^ u{x,t)\. 

Proof. We take the linearized operator L by (|9.4p . as in the proof to Theo- 
rem l9.2l and v{x,t) = G{x, Vu{x,t)) with G{x,p) constructed in Corollary lA. 21 
By the same calculation as ()9.5|) . with G{x,p) in place of F{x,p) we see that, 

Lv < C(l + tr {w'^)) + Gp,p,Wki + Gp.Dp^Bwik 
<C{l + tT{w'')+triw,j)). 



Now since Wij is always positive definite, and det (wij) = -Be" has an upper 
bound by Theorem 18.11 we can see that 



(tr(z«,,))^ - (51 ^fc)^ = (det(«;„))^(^. 



1 



Ai . . . Xk-i^k+i • • • A„ 
<CA^i 
<CtrK^') 

where < Ai < • • ■ < A„ are the eigenvalues of Wij. Since -^^^-rj < C as in the 
proof of Theorem 19.21 we have 

Lv<CtT K^) (l + ^4^1 < Cti- K-')(l + A?)^ 

As in the proof of Theorem [121 define 8 by dHU), only with Ci = C(l + Af)^, 
and note that w(-, t) = for all x G 917 so every such point is a minimum of v. 
Then using the same comparison argument, only in the direction of /3 (which is 
permissible by the obliqueness condition), we obtain 

unp- Afifi=Wfisi^vsi<C{\ + M)'^ (11.1) 

giving the desired estimate. D 
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Theorem 11.2. For any x G dO., and for each t G [0, tmax), 

sup \V^u{x,t)\ < C 
xedn 

for some C that depends on Q, fl* , B, c, and uq, but is independent oft. 

Proof. Assume that ii < tmax and sup sup w^^{x, t) occurs at f = ei 

aox[o,ti] 1^1=1, ({,iy)=o 

and some (a;o,io): and wii(a;o,io) > 1- We write ei — t + h(i where h = -j^x- 
and T — ei — 6/3. Note that (r, u) = 0. Then we obtain, at any x G dO, and any 
i, using Theorem 1 11.11 Theorem 19.21 and the fact that upr ~ 0, 

Wll = Wtt + 2bWrl3 + h'^wpp 

< \T\'^wn{xo,to) + 2bAr0 + b^C{l + M)^ 

< (1 - 2b{l3, ei) + b^\/3\^)wn{xo,to) + 2bArp + b'Cil + M)'-^ 

< (1 - 2b{l3, ei) + C{iy, eif)wn{xo, to) + 2bA,p + C{v, eif{l + M)^ . 

With G constructed in Corollary IA.2I and h constructed in Theorem lA.ll we 
can extend /? and v to all of fl using the formulas i^^x) — \Jh{x) and /3{x,t) = 
Gp^{x,'Vu{x,t)), which also extends b and r to all of il. Since (i'(xo),ei}^ — 0, 
we have V {u{x), ei)'^\x=xo = and thus by Taylor expanding the last term on 
the right hand side above about xq, we have that 

l + 25(/3,ei)-^^^%-<C(l + M)^|x-a;o|2 (11.2) 



Wii{xQ,to) ' Wii{xo,to) 

for all X near xq. 

We now follow a barrier construction as in the proof to Theorem l9.2l Again, 
let fte{xo) = fl n B^{xo) with e chosen small enough so that h < on n^{xo). 

This time we consider the function w(a;,i) = "'^z A —l + 2b(B, ei) !-^4-t + 

KG{x,Vu{x,t)) and 

e(a;,i) = Cih- Ci\x - xqP + v{x,t) - (a, V.yc(a;,yo) -po) 



where po and a are determined the same way as in Tlieorem l9.2l 

First, note that from the positivity of w*^, and writing w as the matrix square 

root of w*^ , we have 

w'^''W^''DiWijDiWik = w''[{DiWij)Wjjn][{DiWik)wkm] > 

Additionally, if we let Mw, M^ be M and M with w in place of u, we find from 
Thcorcm flO.ll that for x e^ and t <ti, 

tr(wy) < CM-u, < C(l + sup|V^u(a;,t)|) 

<C(1+ sup\V^u{x,t)\)^C{l+ sup \wij + Aij\) 
xedn xedn 

< C(l + wii{xo,to)) < Cwii{xo,to) 
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Combining this with (|10.3p . we have that 



'wii{Xo,to) 



-^ZT^T^T ^ -C'(l + tr K^) + tr (w^ )) 



Noting that 26(/3, ei) — — — ,^'^^ ■. is just a function of the form F{x, Vm), 
apphes to give us 

L (2b{f3, ei) - ^^^^) > -C(l + tr K^) + tr (w,,)). 

\ Wii[Xo,to) J 

Thus by fixing a k large enough, and also using the calculation from (19.61) . we 
find that 

Lv> -C(l+tr(w'J)). 

Now since h and G are nonpositive, and by the choice of the linear term, along 
with pT^ . we find that on d{n^{xo)) 

e <C{l + M)'^\x-xof. 

Finally, choosing Ci large enough, we will find that LQ > while < 
on i9(f7e(a;o)) x {0,tmax] and ^e{xo) x {0}, thus by applying the comparison 
principle and differentiating in the direction of — /3 we obtain (V8, —/3) < or 

-wii/3 < C{l + M)^Wii{xo,to) 

<C{l + M)^^ <C{l + M^)^^. (11.3) 

Now by differentiating the condition G{x,Vu{x,t)) = twice and using the 
strict positivity of Gp|.p^ along with (|11.3p . at {xo,to) we get 



^2^^Wkl\^ < Gp^p^UklUli 



<C(1 + M^,)^^ 



Since 



we see that 



M^ = {wii{xo,to)y < ^(u;fei(a;o,io))^ 

k 

<^{C+\uki{xo,tQ)\f 

k 

< C{1 + M„) + ^Ki(xo,io))' 
fe 
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or 

mJ^ - CM^ - C < 

and thus 

M <C + M^<C 

as desired. D 



12 Long time convergence of solutions to the 
flow equation 

By the uniform C^ estimates on u, we find that our equation is uniformly 
parabolic, and the theory of 3, Chapter 14] gives us C^+" estimates on u, and 
hence a standard argument using the Arzela-Ascoli theorem gives existence of 
a smooth solution u for all times i > 0. 

Now fix some positive to a-nd write v{x, t) — u{x, t) — u(x, t + to), and write 
F{x,p,r) =logdet(ry — Aij{x,p)) — B{x,p). Since 

V = F{x, Vu{x, t), V^ u{x, t)) - F{x, Vu{x, f), V^ u{x, t + to)) 

+ F{x, Vu{x, t), V^ u{x, t + to)) - F{x, Vu{x, t + to), V^ u{x, t + to)) 

we can use the mean value theorem to sec that v == a^-'vij + b'''Vi for some 
functions a'-' and 6*. We calculate that 

»i 

[\/r^^F{x,sVu{x,t) + (1 - s)Vu{x,t + to),sV^u{x,t) 



+ (1 - s)V'^u{x,t + to))] ds 
1 
[s V^ u{x, t) + {l- s) V^ u{x, t + to)~ A{x, sVu{x, t) 



+ (1 - s)\/u{x, t + to))Y^ds. (12.1) 

Now, the equation that Wij satisfies combined with bounds on ii and B give us 
a lower bound on tvwij. Combined with the uniform upper bound on Wij, we 
obtain a strictly positive lower bound on the smallest eigenvalue of Wij , uniform 
in t and x. Since u G C'^^'^{fl x M), we have that 

\\u{-,t)-u{-,t + to)\\c-(n) <Ct^, 

and thus by taking to sufficiently small we can ensure that the convex combina- 
tion (1 — s)u{-,t + to) + su{-,t) is close to u{-, t) in C^{i^) norm. By the uniform 
positive lower bound on the eigenvalues of Wij , this ensures that the integrand 
in (|12.ip remains positive definite, hence the equation for v is parabolic. 
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Additionally, we see that for x 6 dO., v satisfies 
= G{x, Vm(x, t)) - Gix, yu{x, t + to)) 

Gpf, {x, sVu{x, t) + (1 — s)Vu{x, t + to))ds J Vk 



= ■■ a^Vk 

By Theorem 19.21 we have that Gp^^{x,'^u{x,t))vk > C > for some C uniform 
in t and x. Thus as above, by choosing to sufficiently small, we can ensure that 
o''vk> % > ^ and we see that v satisfies a linear, uniformly oblique boundary 
condition. 

Now following j6l Section 6.2] we obtain a translating solution of the same 
regularity as u, ie. a function u°°(a;,t) — u°°(a;,0) + Coo • t for some constant 
Coo that satisfies the equation (|3.2a|) . such that ||u — u°°\\ck — > as t — !■ cx) for 
any 1 < k < A. Thus we find 

eC~/(x) = lim e"'"/(a;) = lim det (V^ u°° - A(j;, Vu°°))B~i(a:, Vm°°) 

i— >oo t—^oo 

= lim det (V^ u - A{x, Vu))B-^[x, Wu) 

t—¥00 

= \im det [DT{x,Vu)]g{T{x, t)). 

f— >oo 

We integrate both sides over f2, the C^ estimates on u along with bounds on the 
derivatives of c, / and g allow interchange of the integral and limit. Using the 
change of variables formula with the mass balance condition (12.11) . we obtain 
that Coo = 0. Hence, u°°{x,t) = u°°(a;, 0) is independent of t, and by the 
convergence of u to u°° in the appropriate C*' norms, we see it satisfies the 
desired elliptic equation p.4ap . while Corollarv 15.51 gives the desired mapping 
condition (j3.4bp . 

This completes the proof of the main theorem 13.11 



Appendices 



A Construction of h* and h 

Here we will show the construction of h* and h necessary to carry out the 
barrier arguments in the body of the paper. These constructions appear to be 
common knowledge in the field, but have not been explicitly written down to 
the knowledge of the author (c.f. (9l Section 2]). 
Define the sets 

Te = {x e rj|dist(a;,arj) <e} 

and 

r*. ={yen*\distiy,dn*) <£*} 
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Theorem A.l. Assume that J7 and J7* are uniformly c and c*-convex with 
respect to each other. Then there exist C^ functions h on Vl, and h* on fl* , and 
constants e, e* , 5q, Sq > satisfying the following properties: 

1. V/i = V on do. 

2. h<Q onVe 

3. [D,Jh{x)-c^^^c^J.l{x,y)Dkh{x)]i,^J > 5o\^\^, Vx £ T^, y £ n* , and ^ e M" 

4. Vh* = V* on dfl* 

5. h* <Q on r^, 

6. [D,,h*{y) ~ c''''cuj{x,y)Dkh*{y)]U3 > ^o*l?^ Vy G T*, , x e Q, and 

CgR" 

Proof We will construct h, the construction for h* is similar, but with the 
variables reversed. 

Fix ay G ft*. Let h{x) ~ Cd^{x) — d{x), where d{x) = dist(x,9ri), and 
C > is a constant to be picked. We calculate that 

D,h ^ {2Cd - l)di 

and 

D,jh = {2Cd - l)d,j + 2Cd,dj. 

Clearly, Vh{x) = -Vd{x) = v{x) for x G dfl. 

Now, fix a point x G dfl, and take any ^ G K". We decompose ^ — t{x) + 
av{x) for some a G M and t{x) tangential to dO. at x. Then, 

[Dijh{x) - c'''''cij^i{x,y)Dkh{x)]£,i^j = [-dij{x) + c^'''cijj{x,y)dkix)]£,i£,j 

+ [2Cd,d,]^,^, 
= 1 + 11. 

Considering the matrix didj, we see it has as a basis of eigenvectors: — Vd = v 
with corresponding eigenvalue |V(i|^ = 1, and n — 1 orthogonal vectors, with 
corresponding eigenvalues of 0. Thus, 

// = 2Cdidj{TiTj + 2Tiavj + a ViVj) 
= + + 2Ca^ 



and also 



/ = [-dij{x) + ch^Cij^i{x, y)dk{x)]{TiTj + 2TiaVj + a^ViVj) 
= [Dii^j{x) - c^'''cij^i{x, y)iyk{x)]{TiTj + 2aTiVj + a^ViVj) 
= III + IV + V. 
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By the assumption (j2.3p of uniform c-convexity, we see that 

III>5i\t\^ 

while 

V > -Coa^ 

for some Cq > depending only on ||c||(73 and Q. Finally, using Cauchy's 
inequality with an s, we obtain 

Combining these, we find for C sufficiently large, whenever x G dQ 

Now, by taking e sufficiently small, we can ensure 1 > 1 — 2Cd > ^ and that 
ft. < in r^ . Using the compactness of 9f2 and the continuity of h and uniform 
continuity of c^'^Cij^i in x and y, and by taking e smaller if necessary, we can 
ensure that for each x G F^, there exists an xq G dVl so that 

J Wi-diji^) +c^'^c^j^i{x,y)dk{x)) - (-^^(xo) + c''''c,j,i(a;o,2/)4(a;o))]C*OI < fel^l' 



16 1 



\ 2C\[ddo{x) - dd3{xo)]Uo\ < niel' 
Combining these, we find for x £ T^ 

[Dijh{x) - c^'^Cijj{x,y)Dkh{x)]^i£,j 

= (1 - 2Cd)[-d,j{x) + c''''c,j,i{x,y)du{x)%i, + pCd.d.Ca:)]^;^^ 

> (1 - 2Cd)[{-d,j{xa) + c^'^c^j^i{xQ,y)dk{xQ)) 

- \{-d.^{x) +c'-'^Cij^i{x,y)dk{x)) - {-d.^{xn) + ch^Cij^i{xt3,y)dk{xo))\\£.iij 
+ 2C[didj{xQ) - \{didj(x) - di(ij(a;o))|]6^j 

> (1 - 2Cd){-d,j{xQ) + c''''cy,;(a;o,2/)4(a;o))CjCj + 2Cddj(xa%ij 

Now we have 

\VI> {~dij{xo) + c^'''cijj{xo,y)dkixo))^iS.j, i-dij{xo) + c'-'''ciju{xo,y)dk{xoj)S,iS,j <0 
\VI > ^{-dij{xo) + c'"''c,jj{xQ,y)dk{xo))£,i^j, (-djj(xo) + c'''''c,jj{xQ,y)dk{xo))£,i£,j > 
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so in either case we obtain 






Si 1 
— > -( 



-M 8'^'^' 

^so\e 



for <5o - ^. 



D 



We now make an auxiliary construction, following [I, Appendix A]. 



Corollary A. 2. There exists a function G{x,p) defined on fix Va;c(f2, fi*) such 
that 

for all ^ G M", for some 62 > 0, depending only on f2* and c. 

Also, for each t £ [0,tmax), G{x,Vu{x,t)) = h*{Y{x,'\/u{x,t))) for all x in 

some neighborhood of dQ (possibly depending on t). 



Proof. Let (j> he a, C^ function on M such that 



0">O 



s > 



w 



Note that we may choose (p so that \(p'\ < 1 everywhere. Then define 

^/ fe'(y(a;,p))-/ii(p) 



Gix.p) 



h'{Y{x,p)) + hi{p) 
2 



Mip) 



Y{x,p)er:, 

else 



where 






for some Ci > to be determined, and K > sup|Vu| (which bounded indepen- 
dent of t by Theorem I6.ip . 
Calculating, we find that for Y{x,p) G F*. , 



Go 



+ 



,fh*iY{x,p))-hl{p) 



,(h*{Y{x,p))~hl{p) 



D^,,,h*iY{x,p)) 
DpkPiKip) 



h*{Y{x,p))-hl{p) 



4" V 2 

■{D,,h*{Y{x,p))~D,,hl{p)) 

Iki+IIki+IIhi. 



iDp,h*{Y{x,p))-D,,hUp)) 
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Since (/>" > 0, we have that Illki is positive semi-definite. 
By a simple calculation, 

T * n . k m -J \ U"^ JfTT- , i jn , k 



,m.l 



c"''=C™''(/l* 



^p,nm J 



so by Theorem lA.il 

Dp^p,h*{Y{x,p))^kCi > S*Y,\c''''Ck\'' > 2<5*|CP, some constant S;. 

i 

Now, since |(/)'| < 1, we have that 

{hi+IhOUi > ^Dp,p,h*{Y{x,p))^k^i > 6m' 
or 

{hi+IIki)^k(i > lDp,p,hlip)Ck^i > i-|^|2 

and hence, Gp^p^ is miiformly positive definite for Y{x,p) G F*. . For (x,p) G 
Int{Y^^{r*,-)), G — hi so Gp^pi is clearly also uniformly positive definite there. 
We will conclude the proof of positivity by showing that h*{Y{x,p)) < hl{p) 
on 5(y~i(F*.)) n Int{Ti x V^c{Ti,TF)), hence G = hi on a neighborhood of 
that set. We claim that if (xo,po) e d{Y'^{Tl,)) n Int{n x V:,c(fl, IF)), 
then d(y(a;o,po), (^^*) — £*■ Indeed, we can just use the continuity Y and the 
distance to the boundary, along with taking sequences lying in Y^^(T*t) and 
Int{Q X Vxc{fl, r2*))\F~^(r*. ) that converge to {xo,po). With this, we can see 
that h*{Ylxo,Po)) = (C(e*)2 - e*) = e*{Ce* - 1) < -^ for e* smah enough 
depending only on C. Now by choosing -^ sufficiently small, we can ensure 

that -^ < -^ < hlip) on d (^"^(F*.)) nlnt{nxVxc{n,n*)), giving us the 
desired claim. This also shows that G is C^. 

Finally, for a fixed t, h*{Y{x,\7u{x,t))) = > /iJ(Vu(a:,i)) for x G dn. 
Thus, by the continuity of Vu and the compactness of dQ, we can find a neigh- 
borhood of ar^ on which Gix, Vu{x, t)) = h*{Y{x, Vu{x, t))) = h*{Y{x, Vu{x, t))). 

D 
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